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ABSTRACT 

A  series  of  photoelastic  experiments  was  performed 

on  plate  specimens  containing  semi-elliptical  surface 

flaws  and  subjected  to  biaxial  bending.  Utilizing  the 

three-dimensional  stress  freezing  and  slicing  technique 

o 

and  the  Schroedl-Smith  method  of  analysis,  Mode  I  and 
Mode  II  stress  intensity  factors  were  determined  at  the 
point  of  maximum  flaw  penetration  and  at  the  point  the 
flaw  intersects  the  plate  surface. 

By  varying  the  ratio  of  the  bending  moments  per¬ 
pendicular  and  parallel  to  the  crack  axis,  the  effect  of 
biaxiality  on  the  stress  intensity  factors  was  observed. 
Results  indicated  that  both  Mode  I  and  Mode  II  stress 
intensity  factors  increased  with  a  corresponding  increase 
in  biaxiality  ratio  and  decreased  with  decreasing  biaxiality 
ratio . 
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CHAPTER  I 


INTRODUCTION  TO  PROBLEM 

1 . 1  General  Introduction 

A  significant  amount  of  engineering  research  is 
currently  being  directed  toward  problems  involving  fracture 
of  structural  components.  These  failures  are  very  often 
associated  with  induced  or  inherent  flaws  in  the  material 
which  tend  to  significantly  reduce  the  calculated  strength 
of  the  structure.  By  utilizing  the  stress  intensity  fac¬ 
tor,  K,  a  function  of  applied  load  and  geometry  of  the 
flaw,  it  is  possible  to  predict  the  reduced  strength  of 
the  structural  components  and  therefore  the  entire  struc¬ 
ture  . 

There  are  three  possible  modes  of  deformation 
associated  with  a  crack  as  shown  in  Figure  1.1.  Mode  I 
is  the  opening  mode  where  the  crack  surface  displacement 
is  perpendicular  to  the  loading  plane.  Mode  II  is  known 
as  the  edge  sliding  mode  and  is  characterized  by  the 
surfaces  of  the  crack  sliding  over  each  other  in  a  direc¬ 
tion  perpendicular  to  the  leading  edge  of  the  crack.  Mode 
III  is  the  tearing  mode  in  which  the  crack  surfaces  slide 
over  each  other  in  a  direction  parallel  to  the  leading 
edge  of  the  crack.  A  stress  intensity  factor,  K^.,  K^, 
and  is  associated  with  each  of  these  modes. 


MODE  I,  OPENING  MODE 


MODE  III ,  TEARING  MODE 


Modes  of  crack  deformation. 
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Recent  efforts  have  concentrated  on  the  application 
of  linear  elastic  fracture  mechanics  to  predict  Mode  I, 

II,  and  III  stress  intensity  factors  for  realistic  struc¬ 
tural  configurations  subjected  to  anticipated  loadings. 
Knowledge  of  the  critical  stress  intensity  factor,  K^,  for 
the  materials  assures  reliable  design. 

A  structural  configuration  of  current  interest  in¬ 
volves  determination  of  the  stress  intensity  factors  asso¬ 
ciated  with  semi-elliptical  surface  flaws  in  plates  sub¬ 
jected  to  uniaxial  bending.  The  configuration  and  loading 
realistically  model  conditions  causing  catastrophic  failures 
in  many  engineering  structures  and  especially  those  fabri¬ 
cated  from  high-strength  metal  alloys. 

1 . 2  Purpose  of  the  Investigation 

The  purpose  of  this  investigation  is  to  extend  cur¬ 
rent  design  capability  by  establishing  the  effect  of 
biaxial  plate  bending  and  crack  orientation  on  K^.,  K^, 
and  In  particular,  stress  intensity  factors  were 

determined  for  various  biaxiality  ratios  and  crack  orien¬ 
tations.  A  knowledge  of  the  effect  of  biaxiality  ratio 
and  crack  orientation  on  stress  intensity  factors  will 
greatly  expand  the  practical  use  of  linear  elastic  frac¬ 
ture  mechanics  in  design. 

1 . 3  Scope  of  the  Investigation 

This  study  will  be  broken  down  into  two  primary 
sections.  The  first  series  of  tests  was  performed  using 
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plates  containing  a  centrally  located  semi-elliptical 
surface  crack  oriented  such  that  0=0  degrees,  Figure  1.2. 
The  second  series  of  tests  was  conducted  utilizing  plates 
with  similar  flaws  at  0  =  45  degrees.  Special  consider¬ 
ation  was  given  to  the  selection  of  a  procedure  by  which 
crack  length,  crack  depth,  and  crack  profile  could  be  kept 
nearly  constant,  thereby  reducing  the  number  of  variables 
which  might  affect  the  stress  intensity  factors. 

A  total  of  five  specimens  were  tested  with  the 
parameter  0=0  degrees.  In  each  of  these  tests,  a  posi¬ 
tive  constant  bending  moment,  M  ,  was  applied  to  the 
plate  perpendicular  to  the  line  of  the  crack.  A  second 
bending  moment,  M  ,  parallel  to  the  direction  of  the 
crack  was  varied  from  =  +Mxx  to  =  ~MXX-  The  actual 
test  values  were  0  =  +1.0,  +0.5,  0.0,  -0.5,  and  -1.0,  where 

<J>  =  M  / M  and  directions  of  M  and  M  are  illustrated 
yy  xx  xx  yy 

in  Figure  1.2.  Mode  I  stress  intensity  factor',  K^,  was 
determined  at  two  points  on  the  crack  front;  Point  A,  which 
is  located  at  a  =  0  degrees,  the  point  of  maximum  flaw 
penetration,  and  Point  B  at  a  =  90  degrees,  the  point  at 
which  the  crack  front  intersects  the  plate  surface.  Both 
points  are  illustrated  in  Figure  1.3.  The  remote  stress, 
aox’  Para^lel  to  the  crack  direction  was  also  determined 
in  the  vicinity  of  Point  B. 

Three  tests  were  conducted  with  0  =  45  degrees. 

The  values  of  $  for  these  tests  were  <2>  =  +1.0,  0.0,  and 
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-1.0.  The  resulting  Mode  I,  Mode  II,  and  mixed-mode  stress 
intensity  factors  were  determined. 

1.4  Approach  Used  in  the  Investigation 

The  investigation  was  entirely  experimental  as 
analytical  solutions  to  these  problems  are  not  available 
at  this  time.  A  three-dimensional  photoelastic  analysis 
utilizing  the  stress  freezing  and  slicing  technique  was 
chosen.  This  technique  will  be  explained  in  greater  detail 
in  Chapter  III. 

Three  methods  of  analyzing  the  photoelastic  data 

can  be  used  to  determine  Mode  I  stress  intensity  factors. 

These  methods  are  critically  reviewed  by  Etheridge  and 

Dally. ^  The  decision  was  made  to  use  the  technique  of 

2 

Schroedl  and  Smith  in  this  investigation.  The  accuracy 
of  this  method  is  predicted  to  be  Kj  ±5  percent  if  the 
fringe  loop  radii  are  measured  without  error.  Use  of 
this  method  also  eliminates  errors  involving  the  deter¬ 
mination  of  the  angle  locating  the  extreme  position  of 
the  fringe  loop . 

One  of  the  loading  configurations  for  9  =  45  degrees 
resulted  in  mixed-mode  fringe  patterns  at  the  crack  tip. 
Smith  and  Smith  devised  a  technique  for  reducing  the  data 
from  these  patterns  and  obtaining  both  the  Mode  I  and 
Mode  II  stress  intensity  factors.  Their  method  will  be 
used  in  this  paper.  Another  load  case  produced  pure  Mode  II 


/--vV V 


I 

/ 


deformation.  The  resulting  fringe  loops  were  analyzed  by 
utilizing  the  stress  equations  developed  by  Irwin^  for 
this  mode. 


CHAPTER  II 


REVIEW  OF  RELATED  LITERATURE 

2 . I  Surface  Flaw  in  Tension 

During  the  first  forty  years  following  Griffith's 
findings^  that  the  actual  fracture  stress  of  brittle  ma¬ 
terials  was  lower  than  the  theoretical  values  calculated, 
all  fracture  analyses  were  based  on  through- the -thickness 
cracks.  This  was  due  to  the  difficulties  encountered 
when  attempting  to  analyze  the  surface  crack  which  is 
inherently  three-dimensional.  The  first  contribution  to 
the  problem  of  the  surface  flaw  subjected  to  remote  ex¬ 
tension  was  made  by  Irwin^  who  derived  an  approximate 
expression  for  the  stress  intensity  factor  for  any  point 
on  the  edge  of  a  semi-elliptical  surface  crack.  This 
derivation  was  based  on  Green  and  Sneddon 's^  solution  for 
a  flat  elliptical  crack  in  an  infinite  solid,  Wigglesworth' s 

Q 

solution  of  an  edge-cracked  semi-infinite  solid,  and 
three  corrections .  The  first  correction  accounted  for 
the  addition  of  a  free  surface  normal  to  the  crack  through 
the  major  axis  of  the  ellipse.  The  second  correction  con¬ 
sidered  increasing  the  crack  length  by  an  amount  equal  to 
the  plastic  zone  size  previously  -calculated.  The  final 
correction  was  fcr  the  presence  of  the  free  surface  repre¬ 
senting  the  side  of  the  plate  opposite  the  surface  flaw. 
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The  stress  intensity  factor  at  the  deepest  penetration  of 
a  semi-elliptical  crack  was  given  by: 


where 


Kr  =  l.lo  (^) 


1/2 


(2.1) 


9  _  2 

Q  =  [E(k)]^  -  0.212(-— )  ,  (2.2) 

ys 

E(k)  =  the  complete  elliptical  integral  of  the 
second  kind, 

a  =  the  applied  stress  , 

and 

OyS  =  the  yield  strength  of  the  material. 

This  solution  was  later  found  to  be  inadequate  for  flaws 
which  are  deep  in  comparison  to  plate  thickness  and  frac¬ 
ture  in  the  presence  of  large-scale  yielding.  Early 
attempts  to  estimate  magnification  factors  for  deep  sur¬ 
face  flaws  as  a  correction  to  Irwin's  solution  were  made 

9  10 

by  Kobayashi  and  Smith. 

In  1967,  Smith,  Emery,  and  Kobayashi^  derived  a 

theoretical  solution  for  a  semi-circular  flaw  in  a  half-  • 

space  in  tension.  This  was  extended  to  a  solution  for  a 

part-circular  surface  flaw  in  a  semi-infinite  solid  by 
12  13 

Smith  and  Alavi.  Kobayashi  and  Moss  presented  further 
work  in  the  same  area  including  a  correction  for  the  effect 


11 


of  plasticity.  Thresher  and  Smith, ^  using  the  same 

11  12 

techniques  applied  earlier,  ’  extended  the  work  and 
derived  stress  intensity  factors  for  a  circular  crack 
partially  embedded  in  a  solid  with  finite  thickness  under 
tensile  loading.  Separate  from  the  above  method  of 
approach,  Rice  and  Levy^  obtained  an  approximate  solu¬ 
tion  for  surface  flaws  in  pure  tension  by  replacing  the 
cracked  section  with  a  continuous  spring.  Then  in  1973, 
Shah  and  Kobayashi^  developed  expressions  for  the  stress 
intensity  factor  for  semi-elliptical  surface  flaws  in  a 
finite  thickness  plate  in  tension.  These  expressions  were 
determined  for  various  crack  depth  to  plate  thickness 
ratios . 

2 . 2  Surface  Flaws  Under  Bending 

Expressions  for  the  stress  intensity  factors  for 
surface  flaws  in  bending  fields  are  not  as  numerous  as 
those  for  tensile  fields,  but  nonetheless,  they  have  been 
obtained  by  a  number  of  investigators.  Smith,  Emery,  and 
Kobayashi"^  provided  an  approximate  solution  to  the  pro¬ 
blem  of  semi-circular  crack  in  a  half-space  under  pure 
bending  using  an  alternating  technique.  Smith^  extended 
this  result  to  the  solution  for  a  semi-elliptical  crack 
in  bending  in  a  plate  of  finite  thickness.  Marrs  and 
Smith^  later  conducted  a  series  of  photoelastic  experi¬ 
ments  utilizing  stress  freezing  and  slicing  techniques 


and  compared  their  results  to  Smith's  theory.  It  was 
found  that  experimental  values  ranged  from  12  to  24  per¬ 
cent  higher  at  a  =  0  degrees  and  from  7  to  19  percent 

lower  at  a  =  90  degrees.  Using  the  same  technique  as 
11  12 

Smith,  Smith  and  Alavi  were  able  to  derive  an  expres¬ 
sion  for  the  stress  intensity  factor  for  a  part-circular 
crack.  In  1970,  Rice  and  Levy^  obtained  an  approximate 
solution  using  their  technique  of  a  continuous  line¬ 
spring  model  to  represent  the  bending  resistance  of  the 
surface  flaw. 

18 

Grandt  and  Sinclair  determined  the  stress  inten¬ 
sity  factor  where  the  crack  border  intersects  the  plate 
surface  by  modeling  a  plate  containing  a  semi-elliptical 
surface  crack  as  a  beam  containing  a  spring  in  place  of 
the  crack.  Their  findings,  along  with  analytical  solutions 
by  Smith^  and  Rice  and  Levy,'*’’’  were  compared  to  experi¬ 
mental  data.  Results  indicated  that  the  Grandt- Sinclair 
theory  generally  falls  within  15  percent  of  the  data  for 
a  =  90  degrees.  The  Rice-Levy  theory  was  found  to  under¬ 
estimate  the  stress  intensity  factor  at  a  =  0  degrees  by 
about  15  percent.  Smith's  analysis  for  a  =  0  degrees 
agreed  within  10  percent  of  the  experimental  data  obtained 
by  Grandt  and  Sinclair,  while  results  for  a  =  90  degrees 
were  10  to  30  percent  greater  than  experimentally  deter¬ 
mined  values.  This  is  similar  to  the  findings  of  Marrs 
17  19 

and  Smith.  Shah  and  Kobayashi  obtained  an  approximate 
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solution  to  the  semi-elliptical  surface  flaw  problem  for 
all  values  of  a/ t  by  utilizing  the  solutions  for  a  semi¬ 
elliptical  surface  crack  in  tension  and  for  an  embedded 
elliptical  flaw  approaching  the  free  surface  of  a  semi¬ 
infinite  solid.  This  solution  exhibited  excellent  agree¬ 
ment  with  experimental  results  obtained  photoelastically 

20 

by  Schroedl,  McGowan,  and  Smith.  Using  the  alternating 

method  suggested  by  Thresher  and  Smith, ^  Smith  and 
21 

Sorenson  determined  an  expression  for  semi-elliptical 
surface  flaws  in  bending. 

Schroedl  and  Smith  performed  a  series  of  tests 
using  three-dimensional  photoelastic  techniques  to  obtain 
stress  intensity  factors  near  points  of  maximum  flaw  pene¬ 
tration  and  where  the  crack  border  intersects  the  tensile 
surface  of  the  plate  for  natural  flaws  under  fields  of 
remote  cylindrical  bending.  Results  were  compared  with 
the  analytical  solutions  of  Smith, ^  Rice  and  Levy, ^ 
Grandt  and  Sinclair, ^  and  Shah  and  Kobayashi.^  it  was 
concluded  that: 

1 .  The  Smith  theory  was  in  general  agreement  for 
most  geometries,  but  did  underestimate  for 
a  =  0  degrees. 

2.  The  Rice-Levy  theory  was  not  appropriate  for 
geometries  studied. 


3. 


The  Grandt-Sinclair  theory  yielded  reasonable 
results  at  a  =  90  degrees,  especially  for 
deeper  flaws . 

4.  The  Shah-Kobayashi  theory  gave  excellent 
results  for  a  =  0  degrees. 

2 . 3  Cracks  in  Multiaxial  Stress  Fields 

Until  recently,  most  published  work  dealing  with 
the  stress  analysis  of  cracks  has  dealt  with  flaws  en¬ 
countered  in  a  uniaxial  stress  field.  This  is  a  result  of 
the  complexity  associated  with  analytical  and  experimental 
analyses  of  cracks  in  multiaxial  stress  fields.  Several 
authors  have  approached  the  problem  of  biaxial  stress 

fields  with  regard  to  the  effects  upon  fatigue-crack 

22 

growth  rates.  Kibler  and  Roberts  conducted  experiments 

using  biaxially  loaded  plates  of  aluminum  alloys  6061-T6 

and  6061-T4  containing  centrally  located  through- thickness 

cracks.  They  demonstrated  an  increasing  residual  strength 

with  increasing  lateral  biaxial  tensile  stresses.  Adams' 

23 

study  also  suggested  that  tensile  transverse  stress  can 
decrease  the  cyclic  crack  growth  rate  for  through- thickness 
cracks,  based  on  either  plastic  zone  size  or  crack  opening 
displacement . 

More  recently,  in  a  series  of  fatigue  experiments 
with  center-cracked  plates  of  polymethylmethacrylate  (PMMA) , 
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Radon,  Leevers ,  and  Culver  obtained  results  showing  that 

crack  extension  behavior  is  independent  of  the  stress 

acting  parallel  to  the  crack  plane. 

The  remaining  point  of  view  is  taken  by  Joshi  and 
25 

Shewchuk  on  the  basis  of  the  results  of  tests  on  aluminum 
alloy  plates  containing  surface  flaws .  With  biaxiality 
ratios  of  c^/cr^  =  1.0,  0.86,  and  0.75,  results  indicated 
that  for  a  given  number  of  cycles,  crack  extension  was  the 
greatest  for  the  case  of  and  decreased  propor¬ 

tionally  as  the  ratio  a decreased.  It  should  be  noted 
that  this  is  the  only  investigation  listed  here  which  was 

conducted  using  surface  cracked  specimens.  Oh,  Vardar, 
o  c. 

and  Finnie  D  have  also  shown  that  an  increase  in  degree  of 
biaxial  tension- tension  stresses  results  in  a  faster  fatigue 
crack  growth  rate. 

Several  other  investigators  have  studied  the  effect 

of  biaxial  loading  on  the  Mode  I  stress  intensity  factor. 

27 

Newman  presented  results  for  stress-intensity  correction 

factors  for  through  cracks  emanating  from  a  circular  hole 

in  an  infinite  plate  subjected  to  biaxial  stress.  These 

results  showed  highest  stress  intensity  factors  for  a 

biaxiality  ratio  of  -1.0,  lower  for  0.0  biaxiality  ratio, 

28 

and  lowest  for  a2^a\  =  Grandt  obtained  similar 

results  in  his  studies  of  through- cracked  fastener  holes 
in  biaxial  stress  fields.  Utilizing  the  concept  of  strain 
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intensity  factor,  Hilton  presented  results  which  indi¬ 
cated  that  for  applied  normal  tractions  beyond  the  small 
scale  yielding  range,  the  addition  of  a  component  of  bi¬ 
axial  tensile  loading  decreases  both  the  plastic  strain 
intensity  factor  and  the  plastic  zone  size.  Compressive 

loading  parallel  to  the  crack  direction  had  the  opposite 

30 

effect.  Miller  and  Kfouri  conducted  a  finite  element 
analysis  of  a  square  plate  containing  a  central  through¬ 
thickness  crack  and  subjected  to  biaxial  stresses.  Their 
results  show  the  equi-biaxial  loading  mode  causes  the 
greatest  separation  stress  but  the  smallest  plastic  shear 
ear  and  crack  opening  displacement.  The  shear  loading 
system  induces  the  maximum  size  of  shear  ear  and  crack 
displacement  but  the  smallest  value  of  crack  tip  separation 
stress.  Using  the  relationship  between  crack  opening  dis¬ 
placement  and  stress  intensity  factor,  these  results  would 
indicate  highest  stress  intensity  factors  for  =  -1.0 

and  lowest  values  for  a  biaxiality  ratio  of  +1.0.  Lee  and 
31 

Liebowitz  obtained  results  indicating  the  same  general 

trend  using  a  nonlinear  finite  element  analysis  for  a 

plate  containing  a  through- thickness  center  crack. 

An  analytical  solution  to  the  problem  was  obtained 
32 

by  Smith  using  plane  stress-strain  energy  release  rate. 
His  results,  which  oppose  those  stated  previously,  show 
that  biaxial  tension,  as  compared  to  uniaxial  tension, 
results  in  an  increase  in  strain  energy  release  with  in¬ 
crease  in  degree  of  biaxiality.  Biaxial  compression 
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applied  parallel  to  the  crack  results  in  a  decrease  in 

strain  energy  with  increase  in  biaxiality.  Strain  energy 

release  rate  is  directly  proportional  to  stress  intensity 

factor;  therefore,  these  results  indicate  higher  values  of 

K-j.  for  positive  biaxiality  ratios  and  lower  values  for 

33 

negative  ratios.  Eftis ,  Subramonian,  and  Liebowitz  re¬ 
veal  similar  results  based  on  the  local  elastic  strain 
energy  rate. 


2 . 4  Cracks  Under  Combined  Mode  Loadings 

Most  of  the  experimental  and  theoretical  studies  of 
cracks  have  dealt  with  Mode  I  loadings.  However,  a  small 
number  of  investigations  have  been  conducted  to  determine 
the  effects  of  mixed-mode  loadings. 

o  / 

Wilson  conducted  a  series  of  experiments  with 

center-cracked  plates  subjected  to  a  uniform  tensile  field. 

The  ratio  of  combined  modes  was  controlled  by  the  angular 

orientation  of  the  crack  with  respect  to  the  stress  field 

direction.  Critical  stress  intensity  factors  were  calcu- 

35 

lated  for  various  crack  angles.  Pook  utilized  the  same 
specimen  and  loading  on  plates  with  cracks  at  various 
angles.  Iida  and  Kobayashi  studied  the  crack  propaga¬ 
tion  in  7075-T6  plates  of  similar  configurations  to  the 
two  studies  mentioned  but  under  cyclic  loading.  Stress 
intensity  factors,  and  K^,.  were  determined  by  the 
method  of  direct  stiffness  for  curved  cracks  generated 
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m. 
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from  the  initially  slanted  cracks.  Propagation  rates  in¬ 
creased  in  the  presence  of  loading. 

37 

Shah  used  the  inclined  center-cracked  specimens 
for  combined  Mode  I  -  Mode  II  loading  and  inclined  surface 
flaws  for  combined  Mode  I-II-III  loading.  From  results 
of  the  center-cracked  specimens,  Shah  noticed  that  Mode  II 
has  a  significant  effect  on  the  load  at  which  fracture 
occurs.  Stress  intensity  factors  were  determined  at  the 
surface  of  the  plate  where  Modes  I  and  II  exist  and  at  the 
base  of  the  crack  where  Modes  I  and  III  are  present.  In¬ 
teraction  curves  were  generated  for  critical  values. 

Two  sets  of  photoelastic  experiments  have  been 
reported  in  the  literature  in  which  isochromatic  fringe 
data  have  been  used  to  determine  mixed-mode  stress  intensity 
factors.  In  the  first,  Smith  and  Smith  conducted  stress 
freezing  tests  on  specimens  exhibiting  various  edge  crack 
configurations  in  finite  strips  under  remote  tensile 
loading.  Values  of  K^.  and  were  calculated  using  equa¬ 
tions  developed  for  the  maximum  in-plane  local  shearing 
stress  and  average  fringe  loop  inclination  near  the  crack 
tip.  This  procedure  will  be  explained  in  greater  detail  in 
Section  3.7.3.  Gdoutos  and  Theocaris  developed  a  new 
evaluation  method  of  the  fringe  pattern  for  the  determi¬ 
nation  of  Kj  and  K^.  Extrapolation  laws  were  given  for 
transferring  the  far-from- the-crack- tip  data  to  the  near 


region  in  order  to  bypass  the  error-prone  measurements 
near  the  crack  tip.  This  data  was  then  used  to  solve  the 
equations  noted  by  Smith  and  Smith. 


CHAPTER  III 


EXPERIMENTAL  INVESTIGATION 

3 . 1  Introduction 

The  literature  review  has  revealed  an  experimental 
procedure  for  obtaining  photoelastic  data  for  three- 
dimensional  problems.  It  has  also  illustrated  several 
methods  of  using  this  data  to  calculate  crack-tip  stress 
intensity  factors.  Unfortunately,  the  lack  of  an  analytical 
solution  to  the  problem  under  investigation  does  not  allow 
comparison  between  results  achieved  and  those  that  might 
be  expected.  This  chapter  will  examine  in  detail  the 
models  used,  the  method  of  obtaining  the  data  for  the 
analysis,  and  the  experimental  results. 

3 . 2  Materials  and  Models 

The  material  originally  selected  for  this  investi¬ 
gation  was  PLM-1Z,  obtained  in  liquid  form  from  Photolastic, 
Inc.,  Malvern,  Pennsylvania.  Its  suggested  properties  by 
the  manufacturer  are  given  in  Table  3.1.  It  was  to  be 
cast  in  12  in.  x  12  in.  x  3/4  in.  thick  plates  by  the  in¬ 
vestigator.  However,  after  preparation  of  only  one  plate, 
the  material  became  unavailable  due  to  the  discontinuance 
of  one  of  its  components.  This  necessitated  a  change  to  a 
new  photoelastic  stress  freezing  material,  also  obtained 
from  Photolastic,  Inc.  This  material  requires  substantially 
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Table  3.1.  Properties  of  Photoelastic 
Stress  Freezing  Material, 
PLM-1Z 


Temp .  E  f 

a 

°F  x  1000  psi  psi/fringe/in. 


Coeff.  of 
expansion 
in. /in. /  F 


Poisson' s 
Ratio 


room 
temp . 


400 


80 


7500  39  x  10-6  0.36 


160- 

180  1.8  1.8  >150  90  x  10_D  0.50 


I 

S 
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more  time  to  cast  than  PLH-1Z.  For  this  reason,  it  was 
decided  to  have  the  plates  cast  by  the  supplier.  This 
material  has  not  yet  been  fully  characterized  by  Photo- 
lastic,  Inc.;  therefore,  material  properties  were  not 
given . 

Semi-elliptical  surface  flaws  of  approximately 
the  same  dimensions  as  shown  in  Figure  3.1  were  simulated 
while  casting  the  plates  by  the  use  of  0.005  in.  thick 
shim  stock.  The  end  of  the  shim  stock  was  shaped  to  the 
desired  profile  and  then  suspended  over  the  mold  so  that 
this  shaped  end  produced  an  artificial  flaw.  The  compari¬ 
son  between  this  crack  profile  and  that  of  a  natural 
crack  is  illustrated  in  Figure  3.2.  Model  details  are 
given  in  Table  3.2. 

The  plates  were  machined  to  the  dimensions  shown  in 

Figures  3.3  and  3.4.  Each  plate  had  six  extended  "fingers" 

on  each  side,  excluding  the  zero  biaxiality  specimens 

which  were  only  loaded  on  two  sides.  Plates  were  cut 

using  a  band  saw  with  a  blade  thickness  of  approximately 

0.05  in.  The  dimensions  of  the  remaining  uncut  plate 

were  chosen  as  6  in.  by  6  in.  resulting  in  the  ratio  of 

the  width  of  the  plate  to  crack  length  being  greater  than 
39 

3.0.  Fichter  suggested  that  plates  which  have  ratios 
greater  than  this  value  and  are  subjected  to  a  uniform 
bending  moment  applied  at  its  edges  may  be  treated  as  a 
loaded  infinite  plate. 


x  (in.) 


Figure  3.2.  Comparison  between  natural  crack 
and  semi-ellipse. 


Figure  3.3.  Specimen  configuration  for  4>  =  +1.0,  +0.5,  -0.5,  and 
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3 . 3  Loading  Frame 

The  supporting  body  of  the  loading  frame  was  fabri¬ 
cated  from  1  in.  x  1  in.  x  1/8  in.  steel  angle  sections 
and  1  in.  x  1/8  in.  steel  straps  bolted  together.  The 
actual  loading  mechanism  is  comprised  of  several  parts 
including  the  specimen  support  platform,  scissors  jack, 
and  loading  platform. 

The  specimen  support  platform  was  a  6  in.  x  4  in.  x 
1/4  in.  aluminum  plate  held  by  two  steel  rods  connected 
to  the  loading  platform.  Its  purpose  was  to  support  the 
weight  of  the  model  until  loads  were  applied  at  the  criti¬ 
cal  temperature.  At  that  time,  the  platform  dropped  in 
conjunction  with  the  lowering  of  the  loading  platform, 
allowing  the  model  to  be  loaded  freely.  The  loading  plat¬ 
form  supports  the  weight  which  is  to  be  applied  to  the 
specimen  and  is  lowered  by  means  of  the  scissors  jack. 

Loading  was  achieved  through  the  use  of  40  lb .  test 
mono-filament  fishing  line  and  ball-bearing  fishing  swivels 
The  loading  method  selected  was  similar  to  the  one  used 
by  Rubayi  and  Ved^  and  further  expounded  upon  by  Rubayi.^ 
The  line  was  threaded  through  the  upper  and  lower  frames 
as  well  as  the  model  in  a  manner  illustrated  in  Figure  3.5, 
using  a  continuous  length  of  line  on  each  side  of  the  plate 
The  swivels  were  attached  to  the  bottom  of  the  "fingers" 
by  passing  1  inch  long,  small  diameter  bars  through  the 
eyes.  All  holes  through  which- the  filament  passed  were 
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Connected  to 
Scissors  Jock 


Figure  3.5.  Arrangement  of  loading  the  model 
by  nylon  filament. 
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shamfered  and  corners  rounded  smoothly  in  order  to  reduce 
friction.  To  further  reduce  fricition  and  the  possibility 
of  twist  in  the  line,  all  lengths  were  plucked  and  the 
specimen  was  leveled  by  adjusting  the  ends  of  each  line  by 
the  use  of  small  turnbuckles  fastened  to  the  loading  frame. 

The  calibration  beam  was  suspended  from  the  loading 
frame  and  loaded  prior  to  initiation  of  the  stress  freezing 
cycle.  Loading  was  by  the  four-point  bending  method. 

Beam  dimensions  and  load  application  are  illustrated  in 
Figure  3.6. 

At  this  time,  the  loading  frame  complete  with  the 
unloaded  specimen  and  the  loaded  calibration  beam  were 
placed  in  the  Photolastic  P-2307-M  oven. 

3 . 4  Stress  Freezing  Procedure 

The  stress  freezing  procedure  began  with  the  oven 
at  room  temperature.  It  was  increased  at  a  uniform  rate 
of  5°F/hour  until  the  critical  temperature  of  approximately 
250°F  was  attained.  The  load  was  applied  to  the  model 
which  was  allowed  to  soak  it  this  temperature  for  three 
hours.  The  temperature  was  decreased  at  a  rate  of  2°F/hour 
until  the  oven  temperature  reached  200°F.  From  200°F  to 
140°F,  the  rate  of  temperature  decrease  was  changed  to 
3°F/hour.  The  temperature  was  decreased  from  140°F  at  a 
rate  of  5°F/hour  until  the  oven  temperature  dropped  to 
100°F.  At  this  point,  the  oven  was  shut  off  and  allowed 
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to  cool  to  room  temperature.  Figure  3.7  is  a  photograph 
of  the  loading  frame  complete  with  loaded  specimen. 

3 . 5  Slicing  Plan  and  Technique 

In  order  to  analyze  the  three-dimensional  stress 
state  in  the  model  but  still  employ  two-dimensional 
techniques,  several  slices  were  removed  from  the  speci¬ 
men  to  be  analyzed.  These  slices  consisted  of  a  center 
slice  perpendicular  to  the  crack  and  containing  the 
crack  tip  at  the  point  of  maximum  flaw  penetration  and 
a  surface  slice.  Figure  3.8  shows  the  location  of  these 
slices  with  regard  to  the  model.  Slice  thickness  was 
sufficiently  small  to  allow  the  effect  of  strain  variation 
across  the  thickness  to  be  neglected.  Each  slice  was  cut 
by  a  thin  band  saw  blade  using  air- jet  cooling.  Surfaces 
were  then  sanded  to  a  smooth  plane. 

3 . 6  Photoelastic  Investigation 

3.6.1  Calibration  Specimen.  The  stress  optical  coeffi¬ 
cient,  f  ,  was  determined  for  each  specimen  by  analyzing 
the  calibration  beam  in  a  dark  and  light  field  circular 
polariscope.  The  polariscope  used  was  the  modified 
Columbia  polariscope  which  is  described  in  detail  in 
Appendix  A.  The  calibration  beam  was  placed  in  the  index 
matching  fluid  between  the  partial  mirrors  of  the  fringe 
sharpening  and  multiplication  device  described  in 


3.7.  Loading  frame  with  a  specimen 
under  load. 


Figure 
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Appendix  B.  A  plot  of  fringe  order  versus  distance  from 
the  neutral  axis  was  constructed  to  determine  exact  fringe 
order  at  the  edges  of  the  calibration  beam.  The  value  of 
was  determined  from  the  equation: 


f 

o 


Vl 


(3.1) 


where 

P  =  the  applied  load  to  the  calibration  specimen, 
=  distance  between  loading  points, 
h  =  width  of  the  calibration  specimen, 

=  fringe  order  of  upper  edge , 

and 

=  fringe  order  of  lower  edge . 


3.6.2  Model  Slices.  Slices  of  the  model  specimen  were 
immersed  in  a  bath  of  index  matching  fluid  located  between 
two  partial  mirrors  and  in  the  field  of  a  circular  polari- 
scope.  Fringe  multiplication  was  required  due  to  the  low 
loads  used  to  preserve  small  deformations  in  the  material 
above  critical  temperature.  Figure  3.9  is  a  photograph 
of  the  fringe  sharpening  and  multiplication  unit.  Photo¬ 
graphs  were  taken  of  the  fringe  patterns  for  fringe  sharp¬ 
ening  and  dark  and  light  field  fringe  multiplications .  The 
photographs  were  not  used  for  the  measurement  of  fringe 
radii,  however,  as  a  different  approach  was  employed. 


Figure  3.9.  Fringe  sharpening  and  multiplication  unit 


Figure  3.10.  Telemicroscope  with  fringe  measuring 
eyepiece . 
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Incorporated  in  the  10X  power  eyepiece  was  a  movable  cross¬ 
hair,  the  lateral  travel  of  which  is  measured  by  a  series 
of  graduations  on  the  rotating  knob.  This  device  can  be 
seen  in  Figure  3.10.  Calibration  of  the  movable  cross-hair 
enabled  fringe  measurements  to  be  determined  directly  by 
visual  inspection  of  model  fringe  patterns  and  alleviated 
the  necessity  of  determining  photographic  magnification 
factors.  Data  were  taken  from  the  first  series  of  photo¬ 
graphs,  however,  for  comparison  with  data  obtained  by  the 
movable  cross-hair  method.  Results  suggested  that  the 
degree  of  accuracy  of  the  latter  method  was  equal  to  that 
of  the  photographic  method.  As  a  result  of  these  findings, 
the  new  approach  was  utilized  for  the  remainder  of  the  in¬ 
vestigation.  Typical  photographs  showing  both  unmultiplied 
and  multiplied  fringe  patterns  near  the  crack  tip  for 
several  loading  configurations  are  shown  in  Appendix  D. 

3 . 7  Determination  of  Stress  Intensity  Factors 

3.7.1  Mode  I .  After  careful  review  of  the  three  techniques 
used  for  the  determination  of  Mode  I  stress  intensity  fac¬ 
tors  from  photoelastic  data  (see  Section  1.4),  it  was  de- 
cided  to  use  the  method  of  Schoredl  and  Smith  outlined 
below. 

Beginning  with  the  stress  components  in  a  plane 
normal  to  the  crack  border  in  the  form  for  Mode  I  crack 


extension  suggested  by  Irwin: 
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Smith  simplified  Equation  (3.6)  by  neglecting  a  relative 

O  X 

2n2fo2 

to  - to  obtain: 


Kj  =  (7Tr) 1/2  (2) 1/2  ^  -  a( 


(3.7) 


Kj.  and  aQx  can  be  uncoupled  by  taking  the  difference  be¬ 
tween  Equation  (3.7)  for  the  it^1  and  j  ^  fringe  loops. 
Equation  (3.8)  is  the  result  of  this  operation: 


-i  f  (n.  -n.  ) 
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(3.8) 


K^.  is  computed  from  Equation  (3.8)  for  all  possible 
permutations  of  pairs  of  fringe  loops.  The  average  and 
•standard  deviation  of  these  values  of  Kj  are  determined. 

Kj  average  is  recomputed  after  eliminating  values  of 
which  fall  outside  ±  one  standard  deviation. 

Solving  Equation  (3.5)  for  a  by  use  of  the 
quadratic  formula  gives  a  "positive"  root  of: 
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(3.9) 


Using  the  average  computed  by  the  method  described  pre¬ 
viously,  values  of  a  are  determined  for  each  set  of  (n,r) 
data.  The  same  analysis  performed  in  determining  an 
average  value  is  used  to  determine  an  average  value  for 
aQx.  An  example  of  the  computer  printout  for  this  operation 
is  shown  in  Table  3.3. 


3.7.2  Mode  II.  Following  in  the  manner  of  Irwin,  the 
stress  field  in  a  plane  normal  to  the  crack  border  in  the 
form  for  Mode  II  crack  extension  is : 
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where  the  notation  is  given  in  Figure  3.11.  If  all  measure¬ 
ments  of  r  are  taken  at  9  =  0°,  Equation  (3.10)  reduces  to: 
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Table  3.3.  Sample  Printout  of  K^-  Values 


Crack  Ang 

le 

=  45.0  Degrees 

Biaxiality  Ratio 

=  1.0 

Alpha 

=  0.0  Degrees 

Stress  Fringe  Coef. 

=  2.34 

Slice  Thickness 

.0530 

F(I) 

F(J) 

R(I) 

R(J) 

SIF 

1.50000 

1.33333 

.00830 

.01390 

7.39430 

1.50000 

1.66667 

.00830 

.02040 

9.28035 

1.50000 

1.00000 

.00830 

.03780 

9.48657 

1.50000 

0.83333 

.00830 

.08510 

9.77421 

1.33333 

1.16667 

.01390 

.02040 

12.45822 

1.33333 

1.00000 

.01390 

.03780 

11.04994 

1.33333 

0.83333 

.01390 

.08510 

10.34889 

1.16667 

1.00000 

.02040 

.03780 

9.92777 

1.16667 

0.83333 

.02040 

.08510 

10.32358 

1.00000 

0.83333 

' .03780 

.08510 

10.75225 

1.50000 

1.33333 

.00850 

.01220 

10.28773 

1.50000 

1.16667 

.00850 

.02000 

9.77091 

1.50000 

1.00000 

.00850 

.03970 

9.49522 

1.50000 

0.83333 

.00850 

.09430 

9.72063 

1.33333 

1.16667 

.01220 

.02000 

9.30350 

1.33333 

1.00000 

.01220 

.03970 

9.14305 

1.33333 

0.83333 

.01220 

.09430 

9.54524 

1.16667 

1.00000 

.02000 

.03970 

8.98804 

1.16667 

0.83333 

.02000 

.09430 

9.67088 

1.00000 

0.83333 

.03970 

.09430 

10.48600 

Original 

N 

20 

- 

First  Avg 

;.  SIF 

9.88936 

Standard 

Deviation  = 

1.00898 

New  N 

= 

16 

New  Avg . 

SIF 

9.74599 

Standard 

Deviation  = 

.49799 
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Table  3.3.  (Continued) 


Crack  Angle  =45.0  Degrees 


Biaxiality  Ratio 

1.0 

Alpha 

= 

0.0  Degrees 

Stress  Fringe  Coef.  = 

2.34 

Slice  Thickness  = 

.0530 

F  (I) 

F(J) 

R(I) 

R(J) 

SIF 

1.50000 

1.33333 

.00830 

.01390 

7.39430 

1.50000 

1.16667 

.00830 

.02040 

9.28035 

1.50000 

1.00000 

.00830 

.03780 

9.48657 

1.50000 

0.83333 

.00830 

.08510 

9.77421 

1.33333 

1.16667 

.01390 

.02040 

12.45822 

1.33333 

1.00000 

.01390 

.03780 

11.04994 

1.33333 

0.83333 

.01390 

.08510 

10.94889 

1.16667 

1.00000 

.02040 

.03780 

9.92777 

1.16667 

0.83333 

.02040 

.08510 

10.32358 

1.00000 

0.83333 

'.03780 

.08510 

10.75225 

1.50000 

1.33333 

.00850 

.01220 

10.28773 

1.50000 

1.16667 

.00850 

.02000 

9.77091 

1.50000 

1.00000 

.00850 

.03970 

9.49522 

1.50000 

0.83333 

.00850 

.09430 

9.72063 

1.33333 

1.16667 

.01220 

.02000 

9.30350 

1.33333 

1.00000 

.01220 

.03970 

9.14305 

1.33333 

0.83333 

.01220 

.09430 

9.54524 

1.16667 

1.00000 

.02000 

.03970 

8.98804 

1.16667 

0.83333 

.02000 

.09430 

9.67088 

1.00000 

0.83333 

.03970 

.09430 

10.46600 

Original 

N 

20 

* 

First  Avg.  SIF  = 

9.88936 

Standard 

Deviation  = 

1.00898 

New  N 

= 

16 

New  Avg. 

SIF 

9.74599 

Standard 

Deviation  = 

.49799 

! 
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a  =  0 ,  a  =0 
xx  yy 
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Substitution  of  Equation  (3.11)  into  Equation  (3.3)  yields: 
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max 


KII 
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(3.12) 


which  may  be  combined  with  the  stress-optic  law  to  yield: 
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(3.13) 


Solution  of  Equation  (3.13)  for  K^.^.  gives 


K  _  /9wra/2  nfa 

KII  ~  (2TTr)  Jt~ 


(3.14) 


Values  of  were  computed  for  each  set  of  "n"  and  "r" 
data.  These  values  of  the  stress  intensity  factor  were 
used  in  a  least  squares  regression  analysis  to  determine 
Kjj.  at  the  crack  tip  (i.e.,  r  =  0). 

3.7.3  Mixed-Mode .  For  a  crack  .loaded  in  both  Mode  I  and 
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Mode  II,  Paris  and  Sih  give  the  stress  field  in  the 
vicinity  of  the  crack  tip  as: 
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Substitution  of  these  equations  into  Equation  3.3  yields 
the  following  equation  for  the  maximum  in-plane  shear 
stress : 


T 

max 


- - — T-prUK-rSin©  +  2KTTcos0]2  + 

2(2irr)  1/2  1/\  11 

[K2sin0]2} 


(3.16) 


A  sketch  of  the  isochromatics  in  Figure  3.12 

reveals  that  the  loops  for  the  case  of  mixed- mode  loading 

have  the  same  shape  as  the  Mode  I  fringes,  except  that 

they  are  now  centered  about  a  line  which  is  at  some  angle 

0m  to  the  crack  axis.  For  a  given  value  of  r,  this  angle 

denotes  the  line  along  which  x  has  a  maximum  value.  Smith 

°  m 


3 

and  Smith,  using  a  method  suggested  by  Irwin,  set 
2xm/20  =  0,  which  in  conjunction  with  Equation  (3.16) 
determines  the  angle  0  ,  where  x  has  a  maximum  value. 

The  resulting  equation: 

cos  20m  -  j  =  0  (3.17) 

can  be  solved  for  the  ratio  using  a  value  of  @m 

measured  from  pnotoelastic  experiments .  This  ratio 
coupled  with  Equation  (3.16)  can  be  used  to  determine 
and  K^-j.  for  each  fringe  loop.  By  using  a  least  squares 
linear  regression,  these  values  are  extrapolated  to  deter¬ 
mine  the  stress  intensity  factors  at  r  =  0 . 

3 . 8  Experimental  Results 

The  results  of  this  investigation  can  be  best  pre¬ 
sented  in  two  separate  sections .  Case  I  will  discuss  the 
results  obtained  for  the  five  specimens  with  the  crack 
oriented  at  0  =  0  degrees,  while  Case  II  will  deal  with 
the  three  specimens  having  0  =  45  degrees. 

3.8.1  Case  I.  Table  3.4  presents  a  summary  of  specimen 
details  for  the  five  plates  in  this  category.  A  graphical 
presentation  of  the  values  for  a  =  90  degrees  is  shown 
in  Figure  3.13,  where  normalized  values  of  the  Mode  I 
stress  intensity  factors  are  plotted  versus  the  biaxiality 


Table  3.4.  Slice  Details  for  0=0  Degree  Specimens 


Normalized  Mode  I  stress  intensity  factor  versus  biaxiality 
ratio  for  a  =  90  degrees. 
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ratio,  i.  The  line  drawn  through  the  points  was  determined 
by  a  least  squares  linear  regression  analysis.  The  coeffi¬ 
cient  of  correlation  which  was  calculated  (+0.946)  illus¬ 
trates  the  presence  of  a  very  significant  trend.  It 
suggests  that  a  positive  biaxiality  ratio,  as  compared  to 
uniaxial  bending,  shows  an  increase  in  Mode  I  stress 
intensity  factor  with  increase  in  biaxiality.  Conversely, 
a  negative  biaxiality  ratio  results  in  a  decrease  in  Kj 
with  increase  in  biaxiality.  A  similar  trend  is  recognized 
in  Figure  3.14  for  a  =  0  degrees  with  the  exception  of  the 
unexplained  deviation  at  uniaxial  bending  (i.e.,  $  =  0.0). 
The  fact  that  the  same  trend  is  observed  for  both  a  =  90 
and  0  degrees  solidifies  the  findings.  Theoretical  com¬ 
parison  is  available  only  at  <J>  =  ±0.5  which  is  cylindrical 
bending  for  this  material  because  of  the  value  v  =  0.5  at 
the  critical  temperature.  For  the  slice  containing  a  =  90 
degrees,  comparison  is  made  with  the  theories  presented  by 
F.  W.  Smith^  and  Grandt  and  Sinclair.  ^  Actual  comparison 
of  the  theoretical  values  is  shown  in  Appendix  C.  The 
Grandt-Sinclair  theoretical  point  falls  on  the  linear 
regression  line  and  the  Smith  theory  predicts  a  value 
approximately  14  percent  higher.  This  close  agreement  with 
the  experimentally  determined  value  lends  validity  to  the 
remaining  points .  At  a  =  0  degrees ,  the  Smith  theory 
estimate  of  proved  to  be  20  percent  higher  than  the 
experimentally  determined  curve. 


Figure  3.14.  Normalized  Mode'  I  stress  intentsity  factor  versus  biaxiality 
ratio  for  a  =  0  degrees. 


A  qualitative  check  of  the  accuracy  of  the  loading 
can  be  made  by  a  visual  inspection  of  the  isochromatic 
fringe  pattern.  Eftis  ,  Subramonian,  and  Liebowitzq''5 
illustrated  the  effect  of  various  biaxial  loadings  on  the 
tilt  of  the  characteristic  loops.  A  symmetric  loop  occurs 
only  for  equal  biaxial  tension-tension  loading,  while 
uniaxial  and  negative  biaxial  loadings  cause  a  forward 
tilt  of  the  loops.  For  the  experiments  conducted  here, 
the  $  =  0 . 0  and  -1.0  specimen  results  illustrated  forward 
tilts,  while  for  $  =  +1.0  a  small  backward  tilt  was  observed. 
Therefore,  it  can  be  said  that  the  loading  mechanism  chosen 
did  introduce  relatively  accurate  loads . 

Figure  3.15  is  a  plot  of  the  normalized  remote  stress 
parallel  to  the  crack,  c  ,  versus  the  biaxiality  ratio,  $, 

ua 

for  horizontal  slices.  The  term  "remote  stress"  is  some- 


.15.  Normalized  remote  stress  versus  biaxiality  ratio  for 
a  =  90  degrees . 
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plate  thickness, 

bending  stress,  . 

A  partial  plot  of  this  equation  is  represented  as  a  dashed 
line  in  Figure  3.15.  A  comparison  between  the  theoretical 
values  and  the  measured  experimental  values  reveals  that, 
although  the  same  trend  exists,  experimental  values  conform 
more  to  the  equation: 

6M 

CTox  =  v(l  -  $)  .  (3.19) 

No  explanation  is  given  here  for  this  behavior.  What 
should  be  emphasized  is  that  the  trend  was  reproduced  by 
the  experiments .  No  calculations  were  made  of  for 

vertical  slice,  i.e.,  a  =  0  degrees. 

3.8.2  Case  II.  Three  specimens  had  semi-elliptical  sur¬ 
face  flaws  which  were  oriented  at  45  degrees  to  the  di¬ 
rection  of  the  applied  moments.  Specimen  details  are 
compiled  in  Table  3-. 5. 

The  uniaxial  bending  test  was  conducted  first.  As 
expected,  the  surface  slice  revealed  the  presence  of  mixed¬ 
mode  (I-II)  fringes.  No  fringes  could  be  resolved  for  the 
vertical  slice  even  with  large  orders  of  fringe  multipli¬ 
cation.  A  combination  of  Modes  I  and  II  was  anticipated 


t 

and 


6M 

o 


at  this  location.  The  value  of  0m  (see  Figure  3.12)  was 
measured  and  found  to  be  0^  =  23.5  degrees.  Substitution 
of  this  value  into  Equation  (3.17)  resulted  in  the  equality 

KZ1  =  1.47  KI.  (3.20) 

Values  of  K^.  were  determined  using  Equation  (3.16)  and 
plotted  versus  distance  from  the  crack  tip,  r.  This  is 
illustrated  in  Figure  3.16.  Using  a  least  squares  linear 
regression  analysis,  a  value  of  K-j-  was  determined  at  the 
crack  tip,  i.e.,  r  =  0.  The  resulting  values,  as  shown  in 
Table  3.5  are: 

=  9.990  psi/in. 

=  14.685  psi/in. 

Comparison  is  made  between  these  values  and  that  of  = 

12.085  psi/in.  obtained  from  Specimen  No.  1,  where  0=0 

degrees.  Although  the  value  for  0  =  45  degrees  is 
37 

lower,  Shah  found  that  the  presence  of  has  a  signifi¬ 
cant  effect  on  the  at  which  fracture  occurs . 

Equal  biaxial  tension-tension  loading,  0  =  +1.0, 
resulted  in  pure  Mode  I  isochromatic  fringe  loops.  This 
is  what  was  expected  and  reaffirmed  the  belief  that  the 
loading  frame  was  performing  accurately.  The  computer 
program  utilized  for  all  Case  I  studies  was  implemented 
for  the  determination  of  K^..  Normalized  values  of 
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Kt  =  0.387  at  a  =  0  degrees  and  KT  =  0.619  at  a  =  90 
N  XN  • 

degrees  were  obtained  as  opposed  to  KT  =  0.351  and 

N 

Kt  =  0.797,  respectively  for  the  0=0  degrees  test. 

This  represents  a  decrease  in  K-j.  of  22  percent  for 
the  surface  slice  and  an  increase  in  of  10  percent 
at  the  base  of  the  crack.  Using  the  higher  values  for 
a  =  90  degrees,  it  can  be  noted  that  the  crack  oriented 
at  45  degrees  to  the  loading  would  not  be  as  critical. 

The  final  specimen  tested  was  a  plate  containing 
a  crack  at  45  degrees  to  the  applied  moments  and  loaded 
under  a  negative  biaxiality  ratio,  $  =  1.0.  If  the  load 
was  accurately  applied,  it  was  anticipated  that  the  iso- 
chromatic  fringes  of  the  surface  slice  would  reveal  pure 
Mode  II  loading.  This  is  precisely  what  did  occur.  Using 
Equation  (3.14)  and  the  procedure  described  in  Section 
3.7.2,  a  value  of  =  7.92  psi/in.  was  obtained.  Figure 

3.17  is  a  plot  of  versus  distance  from  the  crack  tip 

and  illustrates  how  this  value  was  determined.  A  much 
larger  was  expected  due  to  the  value  of  obtained 
from  the  $  =  0 . 0  specimen.  There  is  a  strong  possibility 
that  the  method  of  -loading  used  in  this  investigation  does 
not  lend  itself  to  the  acquisition  of  accurate  pure  Mode  II 
stress  intensity  factors.  However,  the  mere  presence  of 
pure  Moae  II  deformation  when  anticipated  alludes  to  the 
accuracy  of  the  method  of  loading.  Isochromatic  fringes 
could  not  be  resolved  in  the  a  =  0  degrees  slice  even  at 
high  orders  of  magnification. 

I 


tress  intensity  factor  versus  distance  from  crack  tip 
1.0  and  a  =  90  degrees. 


CHAPTER  IV 


SUMMARY  AND  CONCLUSIONS 

4 . 1  Summary 

A  series  of  photoelastic  stress-freezing  experiments 
was  conducted  in  an  attempt  to  determine  the  effect  of 
various  biaxial  load  ratios  on  the  crack-tip  stress  in¬ 
tensity  factors .  Plates  containing  semi-elliptical  surface 
flaws  oriented  at  0  and  45  degrees  to  the  applied  load 
were  subjected  to  biaxial  bending  moments  ranging  from 
<J>  =  +1.0  to  <J>  =  -1.0.  Stress  intensity  factors  were  experi¬ 
mentally  determined  at  the  interface  of  the  crack  tip  and 
the  plate  surface,  and  the  point  of  maximum  flaw  penetra¬ 
tion  . 

4 . 2  Conclusions 

Based  upon  the  results  of  the  experimental  analysis 
the  following  conclusions  can  be  drawn: 

1.  For  cracks  oriented  at  0  =  0  degrees,  the  ratio 
of  applied  biaxial  bending  moments  has  a  signi¬ 
ficant  effect  upon  crack-tip  stress  intensity 
factor.  ’  For  both  points  on  the  crack  tip 
studied,  a  positive  biaxiality  ratio,  as  com¬ 
pared  to  uniaxial  bending,  shows  an  increase 
in  Kj  with  increase  in  biaxiality.  Conversely, 
a  negative  biaxiality  ratio  results  in  a  decrease 
in  Kj  with  increase,  in  biaxiality. 
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2.  The  chosen  method  of  loading  is  capable  of 
applying  accurate  biaxial  bending  loads  to 
photoelastic  specimens.  Quantitative  and 
qualitative  comparisons  with  available  theoreti¬ 
cal  solutions  illustrated  good  agreement. 

3.  The  optical  method  of  obtaining  fringe  data 
employed  can  be  used  to  obtain  data  with 
accuracy  comparable  to  photographic  techniques . 

The  photoelastic  stress-freezing  method  has  been 
shown  to  be  a  powerful  tool  in  the  study  of  linear  elastic 
fracture  mechanics.  It  is  hoped  that  work  will  continue 
to  determine  an  accurate  conclusion  to  the  controversial 
subject  addressed  in  this  paper. 
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APPENDIX  A 

CONSTRUCTION  OF  POLARISCOPE 

The  polariscope  used  for  this  investigation 
necessitated  rennovation  due  to  the  addition  of  a  fringe 
sharpening  and  multiplication  unit.  This  device  requires 
a  high  intensity  light  source  which  was  not  initially  pres¬ 
ent  in  the  system. 

The  original  design  consisted  of  a  100-watt  mercury 
vapor  lamp,  two  10-inch  collimating  lenses,  two  quarter- 
wave  plates,  and  two  rotatable  polarizers  arranged  as 
shown  in  Figure  A.l.  The  fringe  sharpening  and  multipli¬ 
cation  unit  was  inserted  between  the  two  quarter -wave 
plates.  This  device  is  explained  in  greater  detail  in 
Appendix  B.  It  was  concluded  thr t  a  250-watt  mercury  vapor 
lamp  would  adequately  fulfill  the  light  intensity  require¬ 
ments.  Through  a  series  of  lenses,  the  light  emitted  from 
the  lamp  is  focused  onto  a  pinhole  which  is  located  at  the 
focal  point  of  the  collimating  lens.  To  incorporate  the 
addition  of  the  larger  lamp  and  the  focusing  lenses,  the 
lamp  housing  was  lengthened.  The  heat  generated  by  the 
250-watt  lamp  necessitated  forced  air  cooling  through  the 
use  of  a  small  blower  directed  at  the  ceramic  lamp  base. 
Figure  A-2  is  a  schematic  representation  of  the  updated 
polariscope . 


Figure  A.l.  Arrangement  of  polariscope  components. 


A  teleraicroscope  with  a  movable  crosshair  (described 
in  Section  3.6.2)  and  a  Nikon  fl.4  lens  was  placed  at  the 
focal  length  of  the  second  10-inch  lens.  Photographs  of 
fringe  patterns  were  taken  through  the  telescope  with  a 
Nikon  F  camera.  Figure  A. 3  is  an  overview  of  the  polari- 
scope . 
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APPENDIX  B 

FRINGE  SHARPENING  AND  FRINGE  MULTIPLICATION 

B . 1  Fringe  Multiplication 

Due  to  the  low  loads  necessary  to  ensure  small  de¬ 
formations,  the  resulting  fringe  pattern  in  the  vicinity 
of  the  crack  tip  did  not  contain  enough  fringes  for 
accurate  evaluation  of  the  stress  intensity  factors.  A 
procedure  known  as  fringe  multiplication  increases  the 
fringe  order  at  every  point  and  produces  an  isochromatic 
pattern  containing  several  times  as  many  fringes  as  the 
ordinary  pattern.  This  effect  is  produced  by  causing  the 
light  to  pass  through  the  model  several  times .  Placing 
partial  mirrors  on  both  sides  of  the  model  yields  multiple 
reflection  of  the  light  and  the  desired  effect. 

The  partial  mirrors  are  inclined  to  each  other  at  a 
small  angle.  A  greatly  exaggerated  view  of  the  effect  this 
has  on  the  light  path  is  shown  in  Figure  B.l.  The  forward 
transmitted  beams  are  converged  by  the  second  collimating 
lens  to  form  a  row  of  images  in  the  focal  plane  at  the 
telemicroscope.  As  can  be  seen  from  Figure  B.l,  the  forward 
transmitted  beams  result  in  odd  multiplication  factors  being 
applied  to  the  isochromatic  fringes. 

It  can  also  be  seen  that  a  certain  degree  of  lateral 
travel  exists  between  the  entry  and  exit  points  of  the  rays. 


> 


74 


The  observed  pattern  is  therefore  no  longer  representative 
of  a  single  point.  Reduction  of  the  lateral  travel  can  be 
achieved  by  rotating  the  entire  mirror  assembly  so  that 
entry  and  exit  points  on  the  two  mirrors  are  opposite  one 
another  for  observation  of  forward  transmitted  beams.  The 
optimum  angle  of  incidence  required  for  this  condition  can 
be  found  by  observing  the  images  formed  on  the  plate  con¬ 
taining  the  pinhole.  This  row  of  even-numbered  images  is 
formed  by  the  first  collimating  lens  converging  the  back¬ 
ward  transmitted  rays.  For  observation  of  a  certain  for¬ 
ward  directed  beam,  the  mirror  assembly  is  rotated  until 
the  source  aperture  is  midway  between  the  images  of  the 
next  lower  and  the  next  higher  even  multiplication  factors 
as  shown  in  Figure  B.2. 

B . 2  Fringe  Sharpening 

If  the  two  partial  mirrors  are  set  parallel  to  the 
model,  the  reflected  and  transmitted  rays  are  colinear  with 
the  incident  ray.  The  intensities  of  the  forward-directed 
beams  are  now  superimposed  resulting  in  the  intensity 
pattern  illustrated  in  Figure  B.2.  This  is  known  as  fringe 
sharpening  because  of  the  resulting  pattern  of  narrow  dark 
and  bright  fringes  on  a  grey  background.  Fringe  sharpening 
results  in  a  more  accurate  determination  of  the  location  of 
isochromatics.  Examples  of  fringe  sharpening  and  fringe 
multiplication  can  be  seen  in  Appendix  D. 


PLATE  CONTAINING  PINHOLE 


Figure  B.2.  Proper  location  of  source  aperature 
lor  optimum  angle  of  incidence  for  a 
multiplication  factor  of  seven. 


Figure  B.3.  Dark  field  intensity  distribution  for  sharpened  isochromatic 
fringes.  After  post  (44). 


APPENDIX  C 


CALCULATION  OF  THEORETICAL  VALUES  OF  Kg 

The  Smith  theory^  and  the  Grandt-Sinclair  theory'*'® 
for  Mode  I  stress  intensity  factors  for  a  semi-elliptical 
surface  flaw  in  bending  both  make  use  of  the  equation: 


KI  = 


Vb 


[ na/Q] 


1/2 


(C.l) 


where 

Mg  =  magnification  factor, 

Og  =  applied  bending  stress, 
a  =  crack  depth, 


and 


$  =  complete  elliptic  integral  of  the  2nd  kind. 

M 

For  the  case  where  the  biaxiality  ratio  =  +0 . 5 

xx 


(cylindrical  bending) : 

a  =  0 . 331  in. , 
Og  =  25 . 808  psi , 


and 


Q  =  1.614. 

Referring  to  graphs  of  Mg  for  various  crack-depth  to  crack- 
length  ratios  found  in  references  10  and  18,  the  following 

values  of  MD  arc  found: 


From  Grandt-Sinclair , 


APPENDIX  D 


ISOCHROMATIC  FRINGE  LOOPS 

Figures  D.l  through  D.3  illustrate  photographs 
of  several  specimen  slices.  Included  are  fringe  sharpening 
and  fringe  multiplication. 


Figure  D.l.  Specimen  No.  1;  $=0.0;  0= 

(a)  a  =  90°,  dark  field,  sharpen,  (b)  a  = 
dark  field.  5  X,  (c)  a  =  90°,  light  field 
(d)  a  =  0°,  dark  field,  sharpen,  (e)  a  = 
dark  field,  5  X,  (f)  a  =  0°,'  light  field, 
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Figure  D.2.  Specimen  No.  5;  $=+1.0;  0  =  45° 

(a)  a  =  90°,  dark  field  sharpen,  (b)  a  =  90°  dark 
field,  7  X,  (c)  a  =  90°,  light  field,  7  X,  (d)  a  ■ 
light  field,  sharpen,  (e)  a  =  0°,  dark  field,  3  X 
(f)  a  =  0°,  light  field,  3  X. 


Figure  D.3.  (a)  ,  (b)  ,  (c)  Specimen  No.  4;  <!>  =  0 . 0  ;  0  =  45 

(a)  a  =  90° ,  dark  field,  sharpen,  (b)  a  =  90°,  dark  field, 
5  X,  (c)  a  =  90°,  light  field,  5  X.  (d),(e),(f)  Specimen 
No.  8 ;  <[>  =  -1.0,  0  =  45°  (d)  a  =  90°,  light  field,  sharpen 
(e)  a  =  90°,  dark  field,  5  X,  (f)  a  =  90  ,  light  field, 

5  X. 
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